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ABSTRACT 

Fundamental aspects of the biochemistry of Type III 
restriction endonucleases remain unresolved 
despite being characterized by numerous research 
groups in the past decades. One such feature is the 
subunit stoichiometry of these hetero-oligomeric 
enzyme complexes, which has important implica- 
tions for the reaction mechanism. In this study, we 
present a series of results obtained by native mass 
spectrometry and size exclusion chromatography 
with multi-angle light scattering consistent with a 
1:2 ratio of Res to Mod subunits in the EcoP15l, 
EcoPI and Pstll complexes as the main holoenzyme 
species and a 1:1 stoichiometry of specific DNA 
(sDNA) binding by EcoP15l and EcoPI. Our data 
are also consistent with a model where ATP hydroly- 
sis activated by recognition site binding leads to 
release of the enzyme from the site, dissociation 
from the substrate via a free DNA end and 
cleavage of the DNA. These results are discussed 
critically in the light of the published literature, 
aiming to resolve controversies and discuss conse- 
quences in terms of the reaction mechanism. 



INTRODUCTION 

Bacteria and archaea protect themselves from potentially 
lethal phage infections by expressing various Restriction 
Modification (RM) systems. These dual-activity enzymes 
inactivate the invading foreign DNA by generating double 



strand breaks, while on the other hand modifying 
(methylating), and therefore protecting the host's 
genome from the endonuclease activity (1). Based on 
their subunit composition, cofactor and recognition site 
requirements and mechanism of action, these enzymatic 
systems can be classified as Types I-IV restriction 
endonucleases (REs) (2,3). The Type III REs are hetero- 
oligomeric assemblies, comprising polypeptides encoded 
by the res (restriction) and mod (modification) genes (4). 
They require two copies of their unmethylated, short 
asymmetric recognition site to be present in the target 
DNA, in an inverted repeat orientation (4,5). DNA 
cleavage is Mg 2+ - and ATP hydrolysis-dependent and 
does not require the presence of the cofactor S-adenosyl 
methionine (AdoMet) [although this can play a role as 
allosteric activator and specificity factor (6)]. The most 
thoroughly characterized Type III RM enzymes are 
EcoPI 51 and the closely related EcoPI both from 
Escherichia coli. Despite several decades of work, there 
are still controversies regarding the mechanism of DNA 
cleavage (7,8). One area of dispute is the subunit stoichi- 
ometry of the endonuclease complex, which we consider 
here. 

The products of the EcoPI 51 and EcoPI mod genes are 
73 and 74kDa proteins, respectively (Table 1), capable of 
recognition site binding (5'-CAGCAG-3' for EcoPI 51, 
5'-AGACC-3' for EcoPI) and N6-adenosine methylation 
using AdoMet as a methyl donor. The Mod subunit when 
expressed on its own is a stable protein, forming a dimer in 
solution and works as an independent methyltransferase 
(9,10). The dimer formation of EcoP15I Mod has been 
recently supported by data from small-angle X-ray 
scattering (SAXS) experiments (11). The K 0 for this 
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Table 1. Theoretical molecular weights (Da) of Res and Mod 
subunits and various subunit stoichiometries for EcoPI 51, EcoPI and 
Pstll 





Res 


Mod 


R,M, 


R,M 2 


R 2 M, 


R 2 M 2 


EcoPI 51 

EcoPI 

Pstll 


110957 
111459 
108264 


74 094 
73 486 
63715 


185051 
184945 
171979 


259 145 
258 431 
235 694 


296 008 
296 404 
280243 


370 102 
369 890 
343 958 



dimerization must be very low (probably in the sub- 
nanomolar range) as free Mod monomers are never seen 
in methyltransferase preparations (12-14). 

The products of the EcoPI 51 and EcoPI res genes are 
nearly identical ~lllkDa proteins (Table 1) which are 
fusions of a Superfamily 2 helicase and a RecB-family 
endonuclease (15,16). A recent study also identified 
several regions within the Res subunit (four in the 
helicase domain and two in the endonuclease domain) 
capable of non-specific DNA binding (17). In complex 
with the Mod subunits, Res forms a holoenzyme 
complex that can hydrolyse ATP and DNA. Attempts to 
express and purify the unmodified Res subunit without 
Mod have proven unsuccessful, most likely because the 
Mod subunit(s) stabilize the Res subunit(s) (18), but a 
fusion protein (with chitin binding protein) has been 
purified (12). 

In an earlier study (12), the molecular mass of untagged 
EcoPI 51 and EcoPI was analysed using analytical ultra- 
centrifugation (AUC) and found to be ~407 and 
~409kDa, respectively. These data were interpreted as a 
Res 2 Mod 2 heterotetrameric stoichiometry (Supplementary 
Figure S3). Quantitative densitometric analysis of samples 
run on SDS-PAGE gels was consistent with this conclu- 
sion (12). A more recent work using EcoP15I enzyme 
supplied by New England Biolabs (NEB) also confirmed 
a heterotetrameric holoenzyme composition using AUC 
and dynamic light scattering experiments (11). In the 
same study the low resolution structure of the EcoPI 51 
endonuclease was determined using SAXS. The data sug- 
gested that the holoenzyme has a crescent shape, most 
likely with the Mod dimer in the centre and one Res 
subunit at each end [similar to the symmetrical arrange- 
ment of HsdR, HsdM and HsdS subunits of the Type I 
REs (19)]. However, as currently no atomic Type III RE 
structure is available, only partial homology models could 
be assembled in the SAXS envelope, indicating a signifi- 
cant level of uncertainty. Although the lowest x 2 values 
were obtained when docking the Res 2 Mod 2 model in the 
envelope, other stoichiometries (e.g. ResjMod^ cannot be 
unambiguously excluded based on these data. 

In contrast to these two studies above, Wyszomirski 
et al. (17) found two oligomeric species in their His- 
tagged EcoPI 51 preparations that could be separated by 
gel filtration chromatography, with Res:Mod ratios of 1:1 
('peak 1') and 1:2 ('peak 2'), as determined by densitomet- 
ric analysis of SDS-PAGE gels. As the latter species 
showed ~6-fold higher specific activity in DNA cleavage 
assays (although DNA binding was not included) and 2- 
to 3-fold higher specific activity in ATPase assays, it was 



considered as the main active EcoPI 51 form. The stoichi- 
ometry of this fraction was studied further by analysing its 
aromatic amino acid composition using second derivative 
ultraviolet (UV) absorption spectroscopy and was found 
to be ResjMod;? (17). This result was also confirmed by gel 
filtration and analytical ultracentrifugation experiments at 
low protein concentrations (<0.4uM), more typical of 
functional Type III RE assays. Furthermore, the AUC 
data suggested that the heterotrimer has an elongated 
shape and/or contains disordered flexible loops. The 
species in peak 1 proved to be inhomogeneous and 
prone to aggregation in the AUC experiments; peak 2 
also showed signs of oligomerization at high concentra- 
tions. It was also reported that when using a low copy 
number plasmid for expression of a non-tagged enzyme, 
the preparations yielded the ReS]Mod 2 form only. The 
authors suggest that the species with 1:1 Res:Mod ratio 
represents a trapped intermediate on the assembly 
pathway of the active heterotrimer. 

Understanding the correct stoichiometry of the Type III 
RM enzymes is important in unravelling their detailed 
mechanism of action. First, there is the question of 
whether one or two Res subunit(s) are required per holo- 
enzyme. In cutting DNA, Type III REs can communicate 
between sites many hundreds of base pairs apart (5,20). 
ATP hydrolysis by the Res subunits plays a key role in this 
communication process. Interpreting how this occurs 
needs clearer knowledge of whether there are one or two 
helicases present. Although earlier models for helicase 
activity suggest a role for dimeric enzymes (21), these 
have been superseded by schemes in which helicases act 
exclusively as either monomers or hexamers/heptamers 
(22). Second, subsequent DNA cleavage requires the inter- 
action of two holoenzymes (12). Thus, there may be either 
two or four endonuclease domains present during the pro- 
duction of a single double-stranded DNA break. 

Third, the subunit stoichiometry is also an important 
question considering the reaction mechanism of methyla- 
tion. Methyltransferases are active as both monomers and 
dimers (23,24). Having two Mod subunits per complex 
raises the possibility of binding two target sequences at 
the same time, giving rise to DNA looping by either the 
isolated methyltransferase or the holoenzyme. Although 
DNA looping is not a necessity for DNA cleavage (25), 
it may play a role in the assembly of a more stable 
complex (26,27). Type III recognition sequences are only 
methylated on one strand (hemimethylation) so dimeriza- 
tion does not play a clear catalytic role in DNA methyla- 
tion. The interaction of two holoenzymes during cleavage 
also brings together four Mod subunits, an apparent 
surfeit given the number of DNA sites (one or two 
depending on the interaction model). However, a 
dimeric methyltransferase may only bind one DNA site 
if either dimerization or DNA binding induces asymmetry 
in the complex. Understanding how the Type III enzymes 
bind DNA would also provide a useful model for 
comparison with other dimeric methyltransferases. 

Due to the conflicting data in the literature on the 
subunit composition of Type III REs, we carried out ex- 
periments using native mass spectrometry (native MS) and 
size-exclusion chromatography (SEC) with multi-angle 
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light scattering (SEC-MALS), techniques that preserve the 
native, untagged complexes and are not sensitive to 
the shape of the molecule. Our data are consistent with 
the heterotrimeric Res 1 Mod 2 stoichiometry. We also show 
that despite two target recognition domains being present 
per enzyme molecule, Type III REs bind only one recog- 
nition site at a time. We discuss our results in light of 
suggested models for the activity of Type III RM enzymes. 

MATERIALS AND METHODS 

Reagents 

Chemicals and reagents were purchased from 
Sigma-Aldrich unless stated otherwise. The term 'apo' 
refers to the nucleotide binding pocket being unoccupied, 
not to the DNA-binding site in the target recognition 
domain. 

Materials, protein preparations and DNA oligonucleotides 

Type III REs were prepared as reported previously (28). 
In brief, wild-type (wt) EcoP15I and wt EcoPI enzymes 
were produced in NovaBlue cells (Novagen) transformed 
with pUC18-based expression vectors coding the entire 
Type III RE operon driven by its natural promoter. 
Following sonication and ultracentrifugation, the crude 
cell extract was passed through a DEAE-Sepharose (ion 
exchange), a heparin (affinity) and finally a Superdex 200 
(gel filtration) column, selecting for fractions enriched in 
the given Type III RE. Preparations typically yielded 10- 
15mg/l of culture. The commercially available EcoPI 51 
was from NEB (catalogue no: R0646S) supplied as a 
10000 U/ml solution in lOmM Tris-HCl pH 7.4, 
100 mM NaCl, ImM DTT, 0.1 mM EDTA, 200ug/ml 
Bovine Serum Albumin (BSA) and 50% (v/v) glycerol. 
The biotinylated EcoPI 51 was generated by engineering 
an Avi-tag to the C-terminus of the Res subunit and co- 
expressing the protein with biotin ligase. The resulting 
labelled protein was affinity purified on SoftLink Soft 
Release Avidin resin (Promega) followed by ion 
exchange using a Mono-Q column (28). The R534A 
mutant of EcoPI 51 was generated by Quickchange muta- 
genesis (Stratagene) using the following oligonucleotides: 
R534A top strand: 5'-CTTCAAGAGGTAGGGGCCGG 
TTTGCGTCTTCC-3' , R534A bottom strand: 5'-GGAA 
GACGCAAACCGGCCCCTACCTCTTGAAG-3' and 
was purified identically to the wild-type. We note that 
all our EcoPI 51 (wt and mutant) and EcoPI preparations 
eluted reproducibly as a single peak from all three 
columns used. Pstll was a kind gift from Alice Sears 
and was expressed in TOP 10 cells and purified on a 
heparin column followed by a Mono-Q column (29). 
SDS-PAGE gels also confirmed that there was no signifi- 
cant amount of degradation products present in our 
enzyme preparations (Supplementary Methods and 
Supplementary Figure SI). Enzyme concentrations were 
calculated using the theoretical extinction coefficients 
according to the Res!Mod 2 stoichiometry. 

Oligonucleotides were purchased from Eurofins MWG, 
with the top strand oligonucleotides either unlabelled or 
5'-labelled with hexachlorofluorescein (Hex) via a 12C 



linker. Two specific DNA (sDNA) duplexes carrying the 
Type III RE recognition sequences (5'-CAGCAG-3' for 
EcoPI 51, 5'-AGACC-3' for EcoPI) and a non-specific 
DNA (nsDNA) duplex were used in this study (Figure 
2 A and Supplementary Figure S4A). The duplexes were 
50 bp in length, carrying the recognition sites six bases 
from the 5' -end where applicable. For annealing the 
duplexes, oligonucleotides were mixed at lOuM concen- 
tration in TE (10 mM Tris-HCl pH 8.0, 1 mM EDTA) at a 
1:1 ratio in the case of unlabelled DNA or at a 1.3:1 in the 
case of labelled DNA. Oligonucleotide mixtures were 
heated to 96°C for 4min and allowed to cool to room 
temperature overnight. 

The assay buffer used throughout the experiments 
(Standard Assay Buffer) was the following: 50 mM Tris- 
HCl pH 8.0, 50 mM KC1, 10 mM MgCl 2 , 1 mM DTT and 
O.lmg/ml BSA. All reactions were carried out at room 
temperature (20 ± 1°C). 

Size-exclusion chromatography with multi-angle 
light scattering 

Protein samples (0.5-1 .0mg/ml, 1.25-2.5 uM) were run in 
Standard Assay Buffer minus BSA at 1 ml/min on a 
Superdex 200 analytical size-exclusion column (GE 
Healthcare), attached to a light scattering diode array 
and a differential refractive index detector (Wyatt 
Technology UK Ltd.). This setup allows for the accurate 
measurement of absolute molecular masses of elutes or 
even a mixture of complexes, independent of shape, 
calculated from the ratio of light scattering and the 
differential refractive index. Chromatograms were 
analysed using the ASTRA software (v6.0.5.3, Wyatt 
Technology UK Ltd.). Where indicated, protein samples 
were incubated with various molar ratios of Hex-labelled 
specific 50-mer DNA duplexes for at least 30min before 
the chromatography run. BSA was omitted from the 
standard assay buffer to avoid possible difficulties with 
peak identification. For samples with 1 mM ATP 
(Sigma), the running buffer was also supplemented with 
the nucleotide. 

To check for reproducibility of the chromatography 
during the SEC-MALS runs, we measured BSA as a 
standard at the beginning and end of each measurement 
day and the elution volumes within the same reaction 
conditions (e.g. in the absence of ATP) were extremely 
reproducible (with an average and standard deviation 
(SD) elution volume of 13.74 ± 0.01 ml). The inclusion 
of ATP in the samples and the running buffer caused a 
slight but reproducible shift (+0.1 ml) in the elution 
volume of BSA. Figures IB, 2B and C and 3 show the 
traces of single chromatography runs rather than 
averages. 

Native MS 

For enzyme-only measurements, samples (10 uM) were 
dialysed against lOOmM ammonium acetate pH 7.0 
prior to analysis. In the case of enzyme-DNA complexes, 
they were mixed in the enzyme storage buffer [10 mM 
Tris-HCl pH 8.0, 100 mM NaCl, ImM DTT and 50 % 
(v/v) glycerol] at a 1:2 ratio (5 uM enzyme and 10 uM 
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Figure 1. Evidence for the heterotrimeric stoichiometry of Type III 
REs. (A) Native mass spectrum of wt EcoP15I. The Y-axis is relative 
intensity, scaled to the most intense peak in the spectrum, which is the 
31+ charge state of the enzyme complex. The numbers above the peaks 
are the measured charge states for the given molecular species. The 
observed molecular weight (259.3 kDa) corresponds closely to that 
expected of the Res^od;, heterotrimer (259145 Da, Table 1). (B) 
SEC-MALS traces of various Type III REs. All enzymes yielded a 
single elution peak with a molecular weight supporting the Res 1 Mod 2 
subunit stoichiometry, except for Pstll, which eluted as a double peak 
indicating a mixture of Mod 2 and Res!Mod 2 species. Solid lines repre- 
sent the normalized light scattering (left Y-axis), whereas dotted lines 
show the calculated molecular masses (right Y-axis). Traces are shown 
in two separate graphs for clarity, noting that the wt EcoP15I and 
EcoP15I obtained from NEB had the same elution volume. See 
"Evidence for a ReslMod2 stoichiometry for Type III REs" in the 
Results section and Table 3 for further details. 



DNA) and magnesium acetate added to the samples at 
2mM. For experiments with ATP, 4mM MgOAc and 
2mM ATP were added. Following 5min incubation, the 
samples were buffer-exchanged into 100 mM ammonium 



acetate using Biorad Micro Bio-Spin P-6 Gel Columns to 
get rid of excess magnesium and other non-volatile buffer 
components that would interfere with the measurement. 
Enzyme stoichiometry measurements were performed on a 
Q-TOF II instrument (Waters) modified for high-mass 
operation (MSVision; including vacuum adaptations, 
32 k quadrupole and extended voltage range) (30,31). 
Commercially available capillaries (PicoTip, New 
Objective) were used for nano-electrospray ionization 
with a spray voltage of 1.6 kV and source temperature 
of 80°C. Gas pressures were 4, 2.3 and 1.2 x 10~ 2 mbar 
for backing, source and collision cell, respectively. 
Tuning parameters were chosen with the aim to preserve 
non-covalent interactions (voltages were 80, 10 and 60 V 
for sample and extractor cone and collision energy, 
respectively). For enzyme-DNA complexes, a Waters 
Synapt G2 HDMS was used with the following param- 
eters: pressures, backing 5mbar, source 5.8 x 10~ mbar, 
trap 4.4 x 10~ 2 mbar; voltages, nano-electrospray using in- 
house-made gold-coated borosilicate capillaries 1.7- 
1.9 kV, sampling and extractor cone 60 V and IV, trap 
and transfer collision energy 50-60 V and 2.1 V, trap DC 
bias 50 V and a source temperature of 30°C. 



RESULTS 

To provide further evidence for the quaternary organiza- 
tion of the Type III holoenzymes, we chose to use two 
techniques not previously applied to these enzymes and 
which can give accurate determination of protein size 
largely independent of overall molecule shape: (a) Native 
MS, which can analyse several species in a single spectrum 
and can provide quite rich detail about complex molecular 
mixtures without the necessity for extensive model fitting or 
assumptions (e.g. as is necessary for AUC) (32-34). One 
shortcoming of this approach is that, for technical 
reasons, the buffer composition is different to the 
standard Type III reaction buffer ('Materials and 
Methods' section). However, it is standard practice in 
native MS to replace other buffers with a compatible, 
volatile buffer such as ammonium acetate, which is chem- 
ically very similar to potassium chloride. The concentration 
and pH of this buffer can be chosen to match the require- 
ments of the protein and other key components (such as, 
e.g. DTT or Mg 2+ ) can be retained. A large body of 
evidence in the literature demonstrates that this MS 
buffer does not typically affect the stoichiometry of 
protein complexes; (b) SEC-MALS equipped with a 
Wyatt light scattering diode array and a differential refract- 
ive index detector, which measures molar mass directly fol- 
lowing separation of different species by conventional SEC 
using Superdex 200 media (35). Although less accurate than 
native MS, SEC-MALS allowed us to use our standard 
reaction buffer, which supports the full enzyme reaction. 
Note that the technique is limited by sample heterogeneity 
following SEC, and the fairly large dilutions during SEC 
which could potentially alter the equilibrium mixtures 
present in the samples. We tested untagged EcoP15I and 
EcoPI samples from a number of sources, including a com- 
mercial stock from New England Biolabs ('Materials and 
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Figure 2. The heterotrimer Type III REs bind a single DNA with a specific site despite the presence of two target recognition domains. (A) Sequence 
of the specific (sDNA) and non-specific ns(DNA) 50-mer DNA duplexes. The 5'-CAGCAG-3' and 5'-AGACC-3' recognition sites for EcoPI 51 and 
EcoPI, respectively, are highlighted with bold and underlined fonts. The cleavage positions (25 and 27 bases downstream of the recognition site in the 
top and bottom strand, respectively) are marked with arrows. (B and C) SEC-MALS traces of a constant amount of EcoP15I (B) or EcoPI (C) with 
sub-stoichiometric amounts to saturating excess of the relevant specific 50-mer DNA duplex (ratios are enzyme:DNA). In the upper graphs, solid 
lines represent the normalized light scattering (left Y-axis), whereas dotted lines show the calculated molecular masses (right Y-axis). In the lower 
graphs, lines show the hexachlorofluorescein (Hex) absorbance traces of the same SEC-MALS runs to confirm the presence of DNA in the eluted 
complexes. (D) Native mass spectra of wt EcoP15I pre-incubated with a 2-fold excess of specific (top panel) and non-specific (bottom panel) 50-mer 
duplex DNA. The Y-axis is relative intensity, scaled to the most intense peak in the spectrum, which is the 9+ charge state of the free DNA. Peaks 
on the left hand side of the dotted vertical line correspond to free excess DNA, whereas peaks on the right hand side can be identified as Res 1 Mod 2 
heterotrimers in complex with a single copy of specific DNA or free Res^od^ holoenzyme in the case of non-specific DNA. As the protein ionizes 
much less efficiently than the DNA, the protein and protein-DNA peaks have been magnified x 1 5 relative to the DNA-alone peaks for presentation 
purposes. 



Methods' section) as well as a biotin-tagged EcoPI 51 used 
in our previous single-molecule studies (28). We also tested 
a sample of Pstll holoenzyme, a more distantly related 
Type III RM enzyme. Molecular weights for different 
theoretical quaternary structures are given in Table 1 for 
reference. To test DNA binding, we utilized 50 bp specific 



DNA duplexes used in previous studies that support 
specific binding and ATP hydrolysis (28) and an equivalent 
non-specific duplex (Figure 2A and Supplementary Figure 
S4A). Theoretical molecular weights for the binding of one 
or two DNAs by different holoenzyme species are given in 
Table 2 for reference. 
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Table 2. Theoretical molecular weights (Da) of various subunit 
stoichiometrics of EcoP15I, the biotinylated EcoP15I and EcoPI in 
1:1 and 1:2 complex with Hex-labelled or unlabelled 50-mer specific 
DNA 



Hex-labelled 50-mer unlabelled 50-mer 







1:1 


1:2 


1:1 


1:2 


EcoP15I 


R,M 2 
R 2 M 2 
M 2 


290658 
401615 
179 701 


322 171 
433 128 
211 214 


289914 
400871 
178 957 


320683 
431640 
209 726 


bio-EcoP15I 


R,M 2 
R 2 M 2 
M 2 


292 989 
406 277 
179 701 


324 502 
437 790 
211214 


292245 
405 533 
178 957 


323 014 
436 302 
209 726 


EcoPI 


R,M 2 
R 2 M 2 
M 2 


289 944 
401 403 
178 485 


321457 
432916 
209 998 


289 200 
400659 
177 741 


319 969 
431428 
208 510 



Evidence for a ResiMod 2 stoichiometry for Type III REs 

To check for the quality of our Type III RE preparations, 
we ran samples on an SDS-PAGE gel, which confirmed 
that the enzymes are highly homogeneous and there is no 
significant proteolytic degradation in the purified samples 
(Supplementary Methods and Supplementary Figure SI). 
Native MS of EcoPI 51, EcoPI and biotinylated EcoPI 51 
showed in all cases almost exclusively one species with a 
Res!Mod 2 stoichiometry (Figure 1A and Supplementary 
Figure S2), in agreement with the enzymes eluting as a 
single peak during purification. For all three enzyme 
complexes a charge state distribution from +27 to +34 
was measured with the most intense peak being always 
the +31 charge state. These charge state distributions 
indicated a similar shape for the different enzymes, as 
the amount of charge a protein accumulates during the 
electrospray process is roughly proportional to the 
solvent accessible surface area for most globular, folded 
proteins (32,36). Experimentally determined masses were 
259.3 kDa for EcoPI 51 and 258.5 kDa for EcoPI, respect- 
ively. The biotinylated EcoPI 51 gave a mass of 261.5 kDa. 
All these values are in close accordance with the theoret- 
ical masses for the ResjMod? heterotrimeric complexes. 

Following SEC of the EcoPI 51 and EcoPI enzymes, we 
reproducibly obtained peaks consistent with single mono- 
disperse species in MALS experiments (Figure IB). The 
determined molecular weights could be most confidently 
interpreted as a ReS[Mod 2 heterotrimer (Table 3). The 
calculated values did not change significantly, and thus 
did not affect the interpretation of the molecular masses, 
when a wider peak selection was used during data 
analysis. Moreover, if Res 2 Mod 2 becomes Res!Mod 2 due 
to sample dilution during the SEC run, there must be 
release of free Res which we would get as a separate 
peak — however we did not observe additional peaks. 
Importantly, our EcoPI 51 preparations of various 
sources showed the same molecular mass as the commer- 
cially available EcoPI 51 from New England Biolabs sug- 
gesting that the fine details of the preparation protocol 
most likely do not affect the stoichiometry of the 



untagged enzyme. We repeated the chromatography runs 
under different buffer conditions (modified NEBuffer4: 
20mM Tris-acetate, pH 7.9, lOmM CaCl 2 , 50 mM potas- 
sium acetate and 1 mM DTT) and got the same results 
(data not shown). For the apo EcoPI 51 in the absence 
of ATP the mean and SD elution volumes for four runs 
were 11.01 ± 0.08 ml, using different enzyme preparations 
including a commercial one from NEB, indicating that the 
chromatography runs were highly reproducible. All the 
data shown are from runs in the Standard Assay Buffer 
to keep reaction conditions consistent across different 
assays. Changes to EcoPI 51, such as the addition of a 
biotin-avidin tag at the C-terminus of the Res subunit 
[used earlier in our single-molecule sliding assay (28)], or 
introducing the R534A point mutation (rendering the 
enzyme ATP hydrolysis deficient) does not influence the 
stoichiometry. We also observed that EcoPI reproducibly 
eluted earlier from the size exclusion column than 
EcoP15I. 

The SEC of the Pstll sample produced two broad peaks 
as observed previously in the absence of Triton X-100 
(37). According to the molar masses derived from the 
MALS, the sample is an equilibrium mixture of 
ResiMod 2 and Mod 2 species, although higher order 
species cannot be ruled out explicitly. The analysis of 
the previous SEC data (37) relied on calibration using 
size standards and on the relative amounts of Res and 
Mod determined using quantitative densitometry follow- 
ing SDS-PAGE. This gave ambiguous results, with a 
retention volume midway between theoretical ReS[Mod 2 
and Res 2 Mod 2 species but with a significant excess of Mod 
over Res. We favoured a Res 2 Mod 2 species based on our 
previous analysis of EcoPI and EcoPI 51 (12), although we 
did not rule out the possibility that the ReS[Mod 2 species 
existed. Our new data does not alter the overall conclu- 
sions of Sears and Szczelkun (37), that subunit assembly 
modulates the activities of Pstll, something that was also 
concluded by Wyszomirski et al. (17) using EcoPI 51. 

Previous sedimentation equilibrium (SE) and sedimen- 
tation velocity (SV) AUC experiments that were inter- 
preted as showing a 1:1 Res:Mod ratio used proteins 
prepared by alternative protocols to those used here 
(11,12). To check how the protein preparations used in 
this study behaved during ultracentrifugation, we revisited 
the SE AUC experiments (Supplementary Methods and 
Supplementary Figure S3). Depending on the protein con- 
centration used and the rotor speed applied, variable mo- 
lecular weights were obtained that are in between that of 
the theoretical ReS[Mod 2 and Res 2 Mod 2 species or, in 
some cases, even above the expected molecular weight 
for the Res 2 Mod 2 species (Table 4). This indicates that, 
probably due to the characteristics of the technique, SE 
AUC gives ambiguous results, possibly reporting on asso- 
ciation/aggregation processes over the long running times 
of the experiment (see 'Discussion' section). SV AUC 
experiments were performed by Janscak et al. (12) previ- 
ously and gave the same 1 : 1 ratio as their SE experiments, 
suggesting that this approach is subject to similar 
limitations. 
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Table 3. Molecular weights (kDa) of free and 50-mer specific DNA complexed Type III REs in the absence 
and presence of ATP as measured by SEC-MALS experiments 
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EcoPI and EcoP15I holoenzyme can only bind one DNA 
molecule at a time 

As all the tested Type III REs contain two Mod subimits, 
the heterotrimeric enzymes could, in theory, simultan- 
eously bind two molecules of DNA containing the recog- 
nition sequence. To look into this, we carried out a 
systematic DNA titration from sub-stoichiometric DNA 
concentrations to saturating DNA concentrations using 
EcoPI 51 and EcoPI and examined the samples first 
using SEC-MALS (Figure 2B and C; Table 3). A Hex- 
labelled DNA was used (Figure 2A) to independently 
monitor the DNA species within the retention volumes. 
ATP was excluded so that only DNA binding would be 
measured (and not motion on DNA or DNA cleavage). 

Following SEC, EcoPI 51 and EcoPI showed different 
profiles as a function of DNA concentration. At sub- 
stoichiometric concentrations, EcoPI 51 seems to be in equi- 
librium between a ResiMod 2 complex bound to one DNA 
molecule (304 kDa) and a less well-defined species with a 
putative molecular mass of ~373 kDa (Figure 2B and Table 
3). This latter species is still visible at equimolar DNA, but 
disappeared as we increased the DNA concentration of the 
specific DNA further. It might represent a state where 
EcoPI 51 undergoes a conformational change upon 
binding to the recognition site and exposes a protein- 
protein interaction surface that can recruit other proteins 
from free solution. However, the estimated molecular 
weight of this species is less reliable due to the sample het- 
erogeneity following SEC and does not correspond clearly 
to a straightforward combination of enzyme(s) and DNA 
based on values in Tables 1 and 2. We suggest that an 
interaction occurs between DNA-bound species and free 
protein(s), although the complexes are clearly heteroge- 
neous based on the light scattering traces. Using excess 
DNA over holoenzyme, only a single species was 
observed with a more defined molecular weight of 306- 
3 1 1 kDa (Figure 2B and Table 3), consistent with one 
molecule of DNA bound to a ResjMod 2 complex (Table 
2). This suggests that, in the absence of free protein, the 
above mentioned interaction does not appear to occur 



between DNA-bound species alone. Under the conditions 
used, we could not drive the formation of a complex with 
two molecules of DNA bound to the heterotrimer. 

In contrast, the different EcoPI complexes formed 
during the DNA titration eluted with the same volume 
(within experimental error) following SEC (Figure 2C). 
A species with a putative average molecular mass of 
~373kDa equivalent to that seen with EcoPI 51 was not 
observed with sub-stoichiometric DNA, suggesting that 
the formation of this complex by EcoPI 51 is most 
probably not an important mechanistic feature of all 
Type III REs. The molar masses derived from the EcoPI 
MALS data are also consistent with only one molecule of 
DNA bound to a heterotrimer (Tables 2 and 3) and the 
Hex traces also confirm the presence of DNA in the 
eluates. We also note that there is a significant amount 
of free DNA present at the sub-stoichiometric DNA 
sample compared with the EcoPI 51 binding experiments. 

In native MS experiments, an enzyme:DNA ratio of 1:2 
was used (Figure 2D). Similar results were obtained when 
the ratio was increased to 1:12 (data not shown). The ex- 
perimentally determined masses were 290.4 kDa for the 
EcoPI 51 DNA complex and 288.1 kDa for the EcoPI 
DNA complex (data not shown). For the biotinylated 
EcoP15I-DNA complex, a mass of 291.0kDa was 
measured (data not shown). These data are all in close 
agreement with the theoretical masses for a heterotrimer 
binding to only a single DNA molecule (Table 2). We note 
that free DNA was visible in the spectra in a larger 
amount than one would have expected at 2-fold excess 
of DNA over enzyme. This may be because the DNA 
ionizes more efficiently than the protein complex and/or 
the binding is less efficient in the native MS buffer than in 
our standard reaction buffer. We also see some free 
Res!Mod 2 and Mod 2 complexes in relatively low 
amounts that are most probably due to degradation of 
the sample during MS. Given that release of a Mod 2 
complex was not observed in the SEC-MALS, we do not 
think that it reflects a normal reaction intermediate. As 
shown in the bottom panel of Figure 2D, the amount of 
free enzyme increased dramatically when non-specific 
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Figure 3. SEC-MALS analysis of the effect of ATP on the protein stoichiometry and DNA binding. SEC-MALS traces of wt EcoP15I (P15), 
EcoP15I helicase mutant (R534A PI 51) or wt EcoPI (PI), with (black lines) or without ATP (grey lines), in the absence (apo) or presence of a 2-fold 
excess of specific 50-mer DNA duplex (1:2). In the upper graphs, solid lines represent the normalized light scattering (left Y-axes), whereas dotted 
lines show the calculated molecular masses (right Y-axes). In the lower graphs, lines show the hexachlorofluorescein (Hex) absorbance traces of the 
same SEC-MALS runs to confirm the presence of DNA in the eluted complexes, where included. The SEC-MALS profiles are not affected signifi- 
cantly by ATP in the absence of specific DNA (A, D). A significant shift in the elution volume and/or a reduction of the calculated mass is indicative 
of a conformational change and dissociation from the specific DNA (B, E). The elution volume of the R534A mutant does not change in the 
presence of ATP and specific DNA (C). Note that for the wt enzymes, the elution volume of the free DNA is significantly increased in the presence of 
ATP (B, E), suggesting that the DNA is cleaved by the enzymes (see "ATP hydrolysis releases the Type III holoenzyme from DNA without affecting 
its stoichiometry" in the Results section and Supplementary Figure S4). 



DNA was used, as only a small amount of enzyme seems 
to bind to this DNA under our experimental conditions. 

ATP hydrolysis releases the Type III holoenzyme from 
DNA without affecting its stoichiometry 

Based on single-molecule and ensemble biochemical 
assays, we have proposed a model for long-range 
interaction by Type 111 restriction enzymes in which 
ATP hydrolysis acts as a molecular switch that causes 
dissociation from the site into a fast DNA sliding mode 
(28). On short DNA oligonucleotides, such as those used 
above (Figure 2A), this rapidly leads to dissociation from 
the DNA by sliding off via the DNA ends. One question 
we sought to address here was, whether the addition of 
ATP would cause the anticipated DNA release and 
whether this process was accompanied by changes in the 
stoichiometry of the holoenzyme. In the absence of DNA, 
wt EcoPI, wt EcoP15I and R534A EcoP15I (mutated in 
Motif VI of the helicase, therefore, lacking ATPase 
activity) had the same heterotrimeric assembly and 
elution volumes following SEC both with or without 
ATP (Figure 3A and D), suggesting that in the absence 
of DNA either these enzymes cannot bind ATP, or that 
binding ATP alone does not affect the subunit stoichiom- 
etry or cause a large conformational change. We note that 
we do not have any direct evidence for ATP binding by the 
wt enzymes in the absence of DNA. Type III REs might 



not bind ATP efficiently in the absence of DNA, as they 
are very weak ATPases when not bound to DNA (5). 

However, in the presence of specific DNA, the elution 
profiles for the wt EcoPI 51 altered upon the addition of 
ATP (Figure 3B), consistent with changes in DNA- 
binding affinity. For EcoPI 51, a broad peak with a 
shoulder was observed with an elution volume intermediate 
of that of the free complex with ATP and the DNA-bound 
complex in the absence of ATP (the shifts being greater 
than the experimental uncertainty in elution volume of 
±0.1 ml — see 'Materials and Methods' section). The popu- 
lation was more heterogeneous and thus the molecular 
mass less accurate, but the calculated range (274- 
280 kDa) is in between that of the free enzyme (264 kDa) 
and the 1:1 enzyme:DNA complex (306 kDa) (Table 3). The 
Hex absorbance profiles show that addition of ATP causes 
release of DNA, although some DNA remains associated 
with the protein peak. This likely represents enzymes that 
have rebound to the site and have yet to undergo the con- 
formational switch. For the R534A EcoPI 51 mutant which 
cannot hydrolyse ATP, no significant change in the elution 
profiles was observed upon addition of ATP (Figure 3C). 
We can therefore interpret the changes seen with wt 
EcoPI 51 as dissociation of the enzymes from DNA follow- 
ing ATP hydrolysis (see below for comments on the change 
in DNA elution volume). We did not observe any peaks 
corresponding to higher or lower order species. 
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Figure 4. Native MS analysis of the effect of ATP on the protein stoichiometry and DNA binding. Native mass spectra of wt EcoP15I pre-incubated 
with specific 50-mer DNA duplex (Figure 2) in the presence of ATP (top panel) or absence of ATP (lower panel). The Y-axis is relative intensity, 
scaled to the most intense peak in the spectrum (the 9+ charge state of the free DNA). Due to lower ionization efficiency, peaks containing protein 
(all peaks rightward of the dotted vertical line) were magnified xl5 for presentation purposes. Fragments containing the intact DNA, cleaved DNA 
and the Res!Mod 2 heterotrimer are labelled with grey, blue and red vertical columns, respectively. In the absence of ATP, most of the enzyme forms 
a stable complex with the specific DNA that partially dissociates to the heterotrimer enzyme and free DNA upon ATP hydrolysis. DNA cleavage 
occurs in the presence of ATP only. The amount of free Res and free Mod 2 increases significantly in the presence of ATP probably due to the 
conformational change associated with sliding (28) that renders the heterotrimer less stable under the native MS conditions. 



For EcoPI, the addition of ATP did not produce a no- 
ticeable shift in the elution profiles (Figure 3D and E). 
Given that the free and DNA-bound species have the 
same elution volume, this is not surprising. However, 
there was a clear loss of the Hex absorbance from the 
protein peak, indicating the release of the DNA. The mo- 
lecular mass calculated from the MALS before and after 
ATP (Table 3) is consistent with release of the DNA 
without any significant change in protein stoichiometry. 
We did not observe any peaks corresponding to higher or 
lower order species. Based on these data, we suggest that 
the Type III enzyme conformation switch following ATP 
hydrolysis is not accompanied by changes in holoenzyme 
stoichiometry. Therefore, the species that dissociates from 
the site and starts to slide is the Res 1 Mod 2 species, as 
expected, based on our single-molecule data where the 
sliding enzyme can rebind the site (presumably because 
the Mod subunits remain associated) (28). 

We observed for both enzymes that the elution volume of 
the free DNA is significantly different to that in the absence 
of ATP (Figure 3B and E, lower panels). One possibility is 
that this is due to the DNA oligoduplex getting cleaved, as 
under the above reaction conditions (high enzyme and DNA 
concentrations and extended incubation times) secondary 
cleavage can occur (6). This was further corroborated by 
the fact that the R534A mutant did not produce such a 
shift (ATP hydrolysis being required for cleavage). To 
examine this further we collected samples following SEC 
and separated the DNA species by native TBE-PAGE 
(Supplementary methods), using both ethidium bromide 
and Hex fluorescence to observe the DNA (Supplementary 
Figure S4B). The results confirm that under these reaction 
conditions, the short single-site DNA is being cleaved by the 



wt enzymes. The kinetics of the cleavage reaction is slower 
than dissociation from the site into the sliding state 
(Supplementary Figure S4C) (also confirmed by Julia 
Toth, unpublished results). We therefore suggest that 
cleavage arises due to an enzyme bound to the DNA 
associating in trans with a recently activated enzyme that is 
still in the sliding conformation (either DNA bound or free). 

To confirm the SEC observations, we carried out native 
MS experiments with EcoPI 51 pre-incubated with 2-fold 
excess of sDNA in the presence (Figure 4 upper panel) or 
absence of ATP (Figure 4 lower panel). In the absence of 
ATP, the majority of the heterotrimeric enzyme formed a 
1:1 complex with a single specific DNA molecule. A small 
amount of free Mod 2 was also observed in the samples as 
well as fragments corresponding to free excess DNA. 
Upon the addition of ATP, characteristic changes were 
noted in the mass spectrum. First, fragments originating 
from free cleaved DNA and the ResiMod 2 heterotrimer 
bound to a single cleaved DNA fragment appeared. 
Moreover, we also observed a marked increase in the 
peaks corresponding to free Res!Mod 2 heterotrimer 
indicative of the partial dissociation of the holoenzyme- 
sDNA complex. We also note that the traces revealed an 
increased level of free Res and free Mod 2 subunits in the 
ATP sample. To check for potential changes in the subunit 
stoichiometry (i.e. the dissociation of the heterotrimer to 
Res and Mod 2 ) following the conformational switch to the 
sliding state we also subjected enzyme samples with sDNA 
and ATP to gel filtration and heparin chromatography 
and ran the resulting fractions on SDS-PAGE gels (data 
not shown). We did not detect any change in the 
Res!Mod 2 heterotrimeric assembly upon ATP hydrolysis. 
In combination with the SEC-MALS data above, we 
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Table 4. Calculated molecular weights (kDa) of wt EcoP15I obtained by analytical ultracentrifugation at multiple protein concentrations at 
various running speeds. The molecular mass values were obtained by fitting a single ideal species model to the experimental data (see also 
Supplementary Methods and Supplementary Figure S3) 

EcoP15I concentration 6 k 7.5 k 10 k Average of 

three speeds 



0.093 mg/ml 349.0 324.7 305.0 318.5 

0.23mg/ml 344.6 325.9 298.5 314.7 

0.46 mg/ml 393.0 392.8 345.6 366.8 

Average of three concentrations 376.0 363.5 330.1 347.6 



suggest that the dissociation of the heterotrimer is a con- 
sequence of the enzyme being less stable in the sliding 
conformation under the native MS conditions. 



DISCUSSION 

As contradicting observations have been published in the 
literature regarding the subunit stoichiometry of Type III 
REs, we carried out native MS and SEC-MALS experi- 
ments to determine the molecular masses of these enzyme 
complexes with improved accuracy. These data strongly 
support a ReS[Mod 2 heterotrimer stoichiometry as 
recently proposed by Wyszomirski et al. (17). We 
did not notice any difference between untagged and 
biotin-tagged protein preparations produced by different 
protocols, labs or strains. We note, however, that tagging 
the protein with an N- terminal His tag and/or changing 
the expression yields by using different promoters can 
produce additional species with subunit stoichiometries 
other than Res!Mod 2 , as seen before by Wyszomirski 
et al. (17). In addition, we were able to determine that 
the holoenzyme only binds a single specific DNA which 
is released following ATP hydrolysis without any change 
in subunit stoichiometry (at least to the majority of 
enzymes under standard reaction conditions). These 
combined results simplify the framework for the future 
consideration of mechanistic models for Type III REs 
(see below). 

The disagreements with previous studies may reflect 
limitations of the experimental techniques used. We 
repeated a sedimentation equilibrium AUC analysis of 
EcoP15I (Supplementary data) and observed significant 
variations in the estimated molecular mass with protein 
concentration and centrifugation conditions (Table 4). 
Similar variations were noted by Reuter et al. and the 
upper end of the estimated values we observed match 
those reported in Janscak et al. (12). We also analysed 
our Coomassie-stained protein gels using densitometric 
analysis. Assuming a uniform dye binding, we applied a 
molecular weight-dependent correction factor and 
obtained results that were consistent with a ResjMod 2 
stoichiometry (data not shown). The simplest explanation 
of the different result seen in Janscak et al. is that the 
earlier protein preparations might have contained 
unbound Res following purification (akin to 'peak 1' 
from Wyszomirski et al.). We would argue that the 
native MS and SEC-MALS data provide a more robust 
analysis as there are a number of assumptions that need to 



be made in analysing the SDS-PAGE gels. The SEC- 
MALS approach described in this article would also be 
appropriate to study the effect of ionic strength, tempera- 
ture and various cofactors (5-adenosyl methionine, S- 
adenosyl homocysteine, sinefungin, etc.) on the subunit 
stoichiometry and/or the ATP hydrolysis-triggered con- 
formational switch to the sliding mode. 

The subunit stoichiometry has important implications 
for the long-range intersite communication process and 
DNA cleavage. Our results clearly demonstrate that 
Type III REs work as monomeric non-classical helicase 
motors. As we did not observe any change in the subunit 
composition following ATP hydrolysis, the sliding compe- 
tent form is the same species as the Res!Mod 2 
heterotrimer capable of recognition site binding. 
Furthermore, according to the ATP hydrolysis switch- 
driven DNA sliding model (28), upon formation of the 
collision complex, two endonuclease domains form the 
active DNA cleavage competent assembly, where a 
single endonuclease is necessary for the cleavage of each 
DNA strand, simplifying possible DNA cleavage schemes 
(12). 

Under the reaction conditions used, Type III REs bind 
only a single copy of a recognition sequence per holoen- 
zyme, despite the presence of two Mod subunits and there- 
fore two Target Recognition Domains. The suggested 
DNA looping states (26,27) must involve two separate 
holoenzymes coming together, one at each recognition 
site, rather than a single enzyme binding two DNA sim- 
ultaneously. Although we did not use a two-site DNA 
substrate in our SEC-MALS and native MS studies, we 
note that we did not obtain any evidence for the formation 
of any higher order species analogous to that in a DNA 
looping complex (e.g. a dimer of ReS[Mod 2 heterotrimers, 
each bound to one DNA). 

Another implication of our observations is that it 
provides information for how dimeric methyltransferases 
may operate (23,38-41). As mentioned above, Type III 
REs bind only one copy of their recognition sites despite 
having two Target Recognition Domains per 
heterotrimeric enzyme, strongly suggesting an asymmetry 
within the two Mod subunits. We did not test the DNA 
binding stoichiometry of the isolated Mod 2 dimer, 
however, we cannot rule out the possibility that the Res 
subunit binding induces the asymmetry. The observed 
DNA binding stoichiometry is also in agreement with 
the enzyme only methylating one strand of the recognition 
site (29,42,43). Other homodimeric methyltransferases 
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can, however, form symmetrical dimers and methylate 
both strands of the DNA in one binding event (as 
opposed to the two events required of a monomeric 
system) (41). The exact nature of the geometry of the 
Type III holoenzyme awaits an atomic resolution crystal 
structure in complex with DNA. 

One of the key aspects of the ATP-driven switch is that 
ATP hydrolysis drives a conformation change into a 
sliding competent state. Binding studies suggest that if 
this structure dissociated from DNA during sliding by 
exiting at a DNA end, the sliding conformation is 
retained for tens of seconds and is unable to rebind 
DNA. We interpret the decreased stability of the 
heterotrimer enzyme in the presence of sDNA upon the 
addition of ATP under the native MS conditions as an 
indication for this conformational difference. 

Finally, in both the SEC-MALS and native MS we 
observed DNA cleavage despite the fact that the rubric 
states that Type III enzymes need two DNA sites in 
inverted repeat on the same DNA for nuclease activity. 
The cleavage seen here does not require direct connection 
between the sites. However, the sliding state must be 
initiated by ATP hydrolysis. We suggest that the sliding 
on DNA acts to not only report on the relative orientation 
of the sites, but elevates the local concentration of 
nuclease domains by a reduction in dimensionality. This 
necessity can be removed by elevating the concentration of 
'activated' enzymes. Under normal conditions, the con- 
centration of enzymes that dissociate from DNA in an 
activated state is too small to produce the protein- 
protein interactions with DNA-bound enzymes. With the 
short single-site DNA duplex used in this study 
(Figure 2A), there is a free DNA end close to the 
cleavage site. We have previously observed that EcoPI in 
an activated form can interact with DNA ends (44). It is 
possible that the elevated cleavage seen here may reflect 
the location of the recognition/cleavage sites. Further 
studies could be informative regarding how activated 
enzymes interact with enzymes bound at recognition 
sites within different substrates. 
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